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ABSTRACT 

Sudden phase anomalies (SPA's) observed on VLF recordings are related to 
Sudden Ionospheric Disturbances (SID's) due to solar flairs. This study presents 
SPA statistics on 500 events observed in New York during the ten year period 
1961 to 1970. Signals were at 10. 2 kHz and 13. 6 kHz emitted from the OMEGA 
transmitters in Hawaii and Trinidad. A relationship between SPA frequency and 
sun spot number was observed. For sun spot number near 85, about one SPA 
per day will be observed somewhere in the world. SPA activity nearly vanishes 
during periods of low sun spot number. During years of high solar activity, 
phase perturbations observed near noon are dominated by SPA effects beyond the 
95th percentile. The SPA's can be represented by a rapid phase run-off which 
is approximately linear in time, peak* tg in about 6 minutes, and followed by a 
linear recovery. Typical duration is 49 minutes. 


INTRODUCTION 

Periodically, portions of the 'surface' of the sun near sunspots are known to 
erupt or 'flare', thereby giving off large bursts of various forms of radiation. 
Flares producing sufficient amounts of energy in the X-ray portion of the spec- 
trum are known to be the cause of sudden changes in the ionization of the earth's 
daytime D-region. This type of Sudden Ionospheric Disturbance (SID) may be 
evidenced as a Sudden Phase Advance or Anomaly (SPA) as observed on record- 
ings of Very Low Freouency (VLF) radio signals received over long paths. The 
typically observed phase variation as a function of time is shown in Figure 1, 
including shojjt term variations due to atmospheric noise; an idealized triangular 
waveform used as a ihfeoretical model for the assumed v; Nation is detailed in 
Figure 2. 


This paper presents the results of a synoptic study of SPA's observed over 
North America from 1961 to 1970 as recorded on VLF signals propagating from 
the Omega Navigation System transmitter in Hawaii to an Omega receiving site 
at Rome, New York. Results of a shorter-term study from 1966 through 1970 
over the path from the Omega transmitter in Trinidad to Rome also are presented. 
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Figure 1. Typical VLF Phase Disturbance Caused by 
Solar X-ray Flare 
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Figure 2. Idealized SPA Shape 
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rates, and cumulative disturbance probabilities. Distribution functions for the 
parameters were obtained together with interrelationships and behaviour with 
carrier frequency, sunspot number, solar zenith angle, and propagation path 
length. 

An Important feature of the present investigation is its restriction to synoptic 
data. No attempt has been made to relate results to particular geophysical mor- 
phologies or to direct data on observed solar X-ray flux. The approach taken 
was direct application to general error considerations for the Omega navigation 
system or other VLF navigation or communication system. For any particular 
specified application, especially one where the X-ray flux from the observed or 
anticipated flares can be determined, more complex analysis may be warranted. 
Additional information on the physics of solar flares especially their effects on 
solar flares can be found in References 1 through 4 of which Reference 4 Includes 
an especially thorough study of one particular SPA for which the associated flux 
data was known. References 5, 6 and 7 provide additional information on the 
Omega system. 

a. Propagation path length 

b. Solar zenith angle on path 

c. Sunspot number 

4. ’Probability’ curves giving the percentage of time that propagation was dis- 
turbed above a certain level. 

Due to the scarcity of Omega data and the low level of solar activity during the 
early 1960's, most of the foregoing results are derived from 10. 2kHz signals 
over the Hawaii to New York path from 1966 to 1970. In general, all results in 
this paper are intended only as summaries of observed effects of SPA's on Omega 
signals and have not been related to possible physical or chemical aspects of 
flare effects. Additional information on the physics of solar flares 1 * 2 * 3 » 4 , as 
well as the Omega System 5 * 6 * 7 , is available in the literature. 


DATA AND COMPILATION PROCEDURES 

The major source of data for this study was the Omega VLF phase difference 
measurements at 10.2 and 13.6 kHz recorded on strip charts at the Rome Air 
Development Center (RADC), 8 New York through 1970. The data represent the 
observed phase difference between the long-path signals from Hawaii (7814 km) 
or Trinidad (3843 km) and the groundwave signal from the local Omega transmit- 
ter at Forestport, New York (36 km). 
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As the signals from all Omega transmitters were being maintained by synchronism 
either by a master-slave mode of operation or by cesium- standard clocks, 9 these 
recordings indicate the absolute effect of ionospheric variation on the behavior of 
the long propagation paths. Data were available essentially continuously from 
January 1966 and paritally from 1961 on 10.2 kHz. 

Statistical data were obtained by manually scanning the recorded charts for SPA's 
and measuring observed times, rates and offsets. The time and phase offset 
determinations were straightforward; rates were indirectly determined by as- 
suming the traingular shape of Figure 2 and obtaining the slope of the onset and 
recovery ’legs’ of the triangle. The onset slope normally followed the recorded 
rapid shift exactly; the recovery slope was determined from an eyeball ’best-fit’ 
line tangent to the half-magnitude point of the recovery period. Deviations from 
the assumed triangular shape were also tabulated for the larger events. Although 
all detectable events were noted, shape characteristics were compiled only for 
those events falling within a 'day window’ centered about mid-path noon. This 
window varied seasonally from 4 hours in winter to 8 hours in summer and was 
typically about 2040 hours per year. This restriction was imposed to prevent 
the curvature in the diurnal phase variation during sunrise and sunset transitions 
from affecting the apparent phase rates-of-change and offsets caused by SPA's. 
Two other restrictions on the compilation of events involved very small (5 cec or 
less) and multiple SPA's. Disturbances of 5 cec* or less have been excluded 
from all analyses due to both the difficulty in detecting such events in the normal 
propagation noise and to the desire to limit the number of events to those exhibi- 
ting a shape sufficient for analysis on the scale of the phase recording charts. 

A 5 cec SPA amounts to less than 0.25 inch on the typical track of an Omega re- 
cording, so detection and/or analysis of smaller events was considered impractical. 

Multiple SPA's; i.e. , additional disturbances occurring before the initial event 
was recovered completely, were included in only those analyses for which the 
relevent shape characteristics could be determined reliably. As an example, 
these events were included in the probability of disturbance and cumulative anal- 
yses but not in any of the shape correlation analyses. The possible effects of 
inclusion or exclusion of multiple SPA's, as well as the other constraints, are 
discussed wherever applicable in the "results" sections to follow. 

Before presenting any results however, the possible effects of equipment and 
the applicability of various types of units will be discussed. 


*1 cec = 1 centicycle = 0.01 cycle. See following sections. 
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EQUIPMENT 


In general, the receiving equipment at the Omega system monitoring sites was 
designed either for use in controlling the timing synchronization of the associated 
transmitters or for general navigation aboard relatively slow-moving vehicles. 

The receiver time constants (approximately 30-60 seconds) and tracking rates 
were, therefore, selected to adequately perform these system-oriented tasks 
with only secondary priority attached to the accurate representation of propaga- 
tion anomalies. Abrupt phase changes were not expected during normal opera- 
tion so that some smoothing over several 10-second transmission format inter- 
vals was not only acceptable but also desirable for noise-reduction purposes. 
None-the-less, the typ’cal time constant was about one minute and the equipment 
was capable of slewing up to lOcec per min, which the following discussion will 
show to be adequate for this study. 

The distortion to input phase variations caused by the receiver can be assessed 
easily by considering two special cases; an abrupt received phase disturbance 
and a ramping one. In the case of an abrupt input phase shift greater than 10 
centicycles (cec), the output would ramp at the maximum slew rate of lOcec/min. 
In practice, a typical magnitude was 15 cec and the typical onset rate was 3.5 cec/ 
minute so that the observed onset rate must be indicative of primarily the actual 
input phase rate. In the case of a ramped input phase variation, the observed 
steady state output will also ramp but delayed from the input by the one-minute 
time constant; within the first minute or so, there will be a slight non-linearity 
betwe n the steady-state ramping and the actual onset. Presumably, this slight 
non-linearity will be difficult to detect due to normal noise effects on the phase 
recordings and the principal effect will be to delay the apparent onset time by 
about one minute rather than affect the observed duration or time to maximum 
onset. Thus the receiving instrumentation is not believed to have had a signifi- 
cant effect on the important statistics considered. 

Some scatter must reflect inability to precisely scale the recorded data. All but 
a few months of data were recorded on linear chart paper having a phase scale 
of 1 cycle full-scale equal to 4.5" and a take-up rate of 3"/hr before 1966 to 
3/4"/hr after 1965. The remaining data were recorded on curvlinear paper with 
a 4.5"/cycle phase scale and a 3/4" /hr take-up rate. The precision of the re- 
corded phase and time readings were about 1 centicycle (cec) and 2 or 3 minutes at 
the slowest chart speed, respectfully. The precision of the slope measurements 
was about 1 degree. As the typical observed onset slope of SPA's was about 85° 
(approx. 3cec/min), the rapid change in the tangent in this region will cause 
the rate uncertainty to increase directly with rate and become approximately 
7 cec /min for 89°. Since the maximum receiver tracking rate of 10 cec /min 
corresponds to angle of 88.5°, the few observed rates in this region and the high 
associated uncertainty lead to the conclusion that th^ equipment was not a factor 
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in limiting the recorded shape behavior of the SPA'S. The two events in 1969 with 
angles of 89° supposedly are not possible, due to equipment limitations, but are 
most probably due to the measurement uncertainty. 


PHASE ANOMALY SIGNIFICANCE AND UNITS 

All primary statistics in this paper are expressed in units of the experimentally 
observed phase changes or rates over the propagation path being investigated. 
While lacking universality, the approach renders the statistics independent of 
assumed ionospheric models. To assist in making the results more meaningful, 
additional scales have been superimposed in which the observed variations are 
associated with "equivalent" changes. 

Experimental measurements were made in angular units of centicycles (1 centi- 
cycle = 1 cec = 0. 01 cycle). It is often common practice to draw an "equivalence" 
between phase and time by associating one cycle of the carrier with one pe- 
riod. ll At 10.2kHz this yields lcec ~lps while at 13.6kHz, lcec ~3/4ps. 

A similar relation between phase and "equivalent" distance when ranging may be 
obtained by associating one cycle with one wavelength. The wave length at 
10.2kHz is 30km so that lcec ~0.3km distance whileatl3.6kHz, lcec ~0.23km 
distance. (In the hyperbolic phase difference mode usually used with Omega, 
these distance equivalences are approximately halved so that lcec -0.15km.) 
Associated rates follow directly: one centicycle/minute corresponds to a frac- 
tional frequency variation of 1.7 x 10~ 8 at 10.2kHz and 1.3 x 10~ 8 at 13.6kHz. 
Velocity equivalences when interpreted for radial ranging are 1 centicycle/minute 
~18 km/hr at 10.2 kHz or 14 km/hr at 13.6 kHz. Another association which is 
often made is that between phase variation and ionospheric height change. In the 
common calculation, the path length is divided by the time equivalent phase varia- 
tion to obtain a velocity variation. The velocity variation is associated with an 
equivalent change in the height of the earth-ionosphere waveguide necessary to 
match the computed velocity variation. The difficulty with this method is that it 
assumes the phase shift is entirely due to an effective height change rather than 
changes in ionospheric electron density gradient or other process. Further, the 
calculation is valid only if the excitation factor remains constant. An additional 
complication, discussed in a subsequent section, is whether the calculation 
should use the entire propagation path length. Despite the limitations, "equiva- 
lent" height variations 7 have been associated with the data based on a partial 
relative velocity variation with height of 3.2 x 10 Am at 10. 2 kHz and 1.9 x 
10 ~*Am at 13.6kHz corresponding to equivalent height changes of 0.12 and 
0.15km per centicycle at 10.2 and 13.6kHz respectively over the Hawaii-New 
York path and twice as much over the Trinidad-New York path. 
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RESULTS 


Maximum Phase Offset Distributions 

Figure 3 gives the number distributions of the maximum phase offsets pro ' aced 
by SPA’s observed at 10.2 kHz on the Hawaii-New York path from 1966 through 
1970. The data was quantized in 5cec intervals, and include all observed events 
greater than 5cec. Multiple SPA's have been included by considering the appar- 
ent additional offset produced by the subsequent events as being individual events. 
This procedure has had the most significant effect in 1967 and 1970 wherein the 
number of observed events has been increased by 20-25%. For other years the 
effects were 10% or less, and no significant effect on the overall distributions 
was found. 

Several features are apparent from Figure 3. There are fewer SPA’s and they 
are smaller during periods of lower solar activity (1966-1967) than near solar 
cycle maximum (1968-1970). The phase change typically peaks between 10-15 cec 
corresponding to ionospheric phase-height changes of 1.2 to 1.8 km. For SPA’s 
greater than 5 cec, the average magnitude was 23 cec corresponding to height 
changes of 2. 8km. It is notable that in most years, fewer SPA's were observed 
with amplitude between 5 and 10 cec than were observed with amplitude between 
10 and 15 cec. The distribution thus does not appear of the type where there is 
a very large probability of small events such as might be expected if SPA's were 
but rare manifestations of a phenomenon characteristic to normal ionospheric 
formation. Instead, the distributions suggest SPA's and the causual solar flares 
are specific identifiable events which tend to have typical magnitudes character- 
istic of the flare formation physics. 

Onset Phase Rate-of-Change Distributions 

The number distributions for the onset rate-of-phase-change for various years 
are shown in Figure 4. The distributions appear similar and essentially normal 
in each year differing primarily in the fewer events observed during lower solar 
activity in 1966. The median onset rate is 3.5cec/minute which corresponds to 
a fractional frequency shift of 6 x 10 8 equivalent to a radial- range rate of 63 km/ 
hr. The associated rate of lowering of the ionospheric phase height is about 0.4 
kilometer per minute. 

The observed onset rates are generally well within instrumental capabilities and 
thus are believed to indicate true ionospheric changes. The two extreme rates in 
1969 may be the result of measurement uncertainty whereas the near absence of 
rates between 4 and 5 and between 6 and 7 cec /min is due to roundoff of the onset 
angle during scaling. 


449 






MAXIMUM PHASE OFFSET IN CEC It CEC * 1 JISI 


0 S 10 

MAXIMUM IONOSPHERIC. PHAiSE HEIGHT CHANGE (km) 


Figure 3. Maximum Phase-Offset Distributions 
Observed During SPA Events on Hawaii to New York 
Path at 10.2 kHz 
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Figure 4. Onset Phase Rate-of-Change Distributions 
Observed During SPA Events on Hawaii to New York 
Path at 10.2 kHz 
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In this and the following 2 sections, multiple SPA's have been excluded because 
of the difficulty in determining the rates-of-change and duration of individual 
events in the disturbance. 

Disturbance Duration Distributions 


Figure 5 shows the number distributions for duration of SPA disturbances obser- 
ved during the same data sample used to show amplitude distribution in the pre- 
vious section. The plots indicate a peaking effect around 40 minutes or so with 
a decrease in frequency with length of disturbance. Except for the minor occur- 
rance of a few longer SPA's in the latter years, there does not appear to be any 
significant variation of disturbance duration with sunspot number, which nearly 
tripled over the time span being analyzed. 

As will be discussed in the following section on decay rates-of-change, the dis- 
turbance duration may have been significantly shortened by the modeling approxi- 
mation used for the 'recovery' leg of the event. For a typical duration of 40 min- 
utes, an error of 10% is significant but probably unimportant in a practical sense. 

Recovery Phase Rate-of-Change Distributions 

Figure 6 shows the number distribution for the decay rates-of-change of phase 
for the same data sample as previously shown. The peak occurrence is at 
0. 5 cec/min or less for each year and the distributions are not quite normal with 
the heaviest weighting occurring at the slower rates. 

As noted in the data action, the decay rate was determined from an approximate 
linear fit to the normally exponential recovery exhibited by the typical SPA. An 
investigation of the possible errors introduced by this approximation indicated 
that the decay slope estimation usually was biased toward the early portion of the 
recovery period and therefore resulted in a higher indicated overall decay rate, 
and subsequently, a shorter indicated duration of the disturbance. This system- 
atic error will not significantly affect the decay rate statistics, as long as they 
are understood to apply to all but the final portion of the recovery period. How- 
ever, the duration statistics may be significantly affected as was discussed 
previously. 

Disturbance Leve l Probability Analysis 

A most useful statistic describing the effects of SPA's is the percentage of time 
that an Omega Phase track is likely to be disturbed above a given value. Such 
information may be extracted from historical data summarized according to time 
span and the percentage of time the observed disturbances exceeded various in- 
cremental levels. Assuming that the chronological variation of SPA occurrence 
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can be related to some solar activity Indicator (e.g. , sunspot number), the his- 
torical statistics may be interpreted as the probability of occurrence of distur- 
bances during any given period of solar activity. 


The historical summaries may be obtained from the data compiled for this re- 
port by assuming SPA's to have the triangular shape of Figure 1 and a reference 
level A <(> c as shown in Figure 7. The total duration T of the SPA is computed 
from the tabulated maximum amplitude A# m and the onset and recovery angles 
6 0 and 0 R : 


T = T 0+ T r 


A* 


m 




Tan d 0 Tan (180 -0 R ) 


* A0 


m 


yTan d Q (Tan 0 R | J 


T = N> 


m 


/Tan 6 0 + |Tan 0 R | \ 
\ Tan 0 O (Tan 0 R I / 


(I) 


Using similar triangles, the duration T c for which A0 exceeds the reference 
level A 4> c can be found from: 
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As written, Equation 2 utilizes the amplitudes and rates of change obtained from 
the disturbances observed over a specific path under specific solar illumination 
conditions. When Equation 2 is divided by the total time available £T 0 for detect- 
ing SPA's, the result is the percentage of time that the observed disturbances 
exceeded the specified reference level: 


T c (%) 



( 3 ) 


If T c (%) is thought of as a percentage probability of occurrance, the notation can 
be changed to read: 


P(A0>A0 c ) 


1 00 T c (%) = 100 


ST: 


2T n 



(4) 


Figure 8 gives plots of P for the years 1966 through 1970 over the range of appli- 
cable reference levels. The relative position of the curves tends to follow the 
variation of solar activity over this period, except for 1967, which appears more 
active than expected. This may have resulted from the procedure used to nor- 
malize the observed disturbed time to the total time available for observation. 
This normalization assumes a random distribution of SPA activity throughout 
the year and neglects any compression of activity into specific months or even 
weeks, which is typical of recent solar activity levels. As 1967 had significantly 
less observed time then other years, any error introduced by the normalization 
would be greater for that year. In addition, 1967 contained a greater number of 
multiple SPA's, which amounted to a significant proportion of the already very 
small percentage of disturbed time. After consideration for the high sensitivity 
of the computation in the region of low percentages, the overall result is in good 
agreement with that expected from the sunspot number variation during this pe- 
riod. The following section discusses the correlation of SPA activity and the 
yearly average sunspot number. 


Frequency of SPA Occurrence with Svn^pot Number 

Figure 9 shows the observed variation of SPA occurrence with the National 
American sunspot number, R A * 2 , The plotted data points have been determined 
from measurements at 10.2 kHz on the Hawaii to New York path and include all 
events greater than 5cec as observed in a typically 6-hour duration 'window' near 
mid-path noon. The points shown for N 24 are an extrapolation to the event levels 
expected under the condition of continuous solar influence over 24 hours. The 
curve for N 12 is a similar extrapolation to a typical 'daylight' condition of 12 
hours somewhere on the globe. 
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Figure 8. Percent Probability of Disturbance Level Occurrence for Observed 
SPA Events on Hawaii to New York Path at 10. 2 kHz 
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Figure 9. Variauon of SPA Occurrence With Sunspot Number 


The effect of multiple SPA's is shown by the double values plotted for certain 
years. The upper point includes all events discernible within the duration of the 
multiple disturbances; the lower point assumes each multiple event to be only a 
single disturbance. The largest effect is in 1967 which exhibited about 20% less 
observation time but proportionally more multiple events than other years. 
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The distribution functions of Figure 3 suggest that perhaps the event levels might 
be about 10% higher, if small comparatively unimportant events of less than 5 cec 
also had been included. 

The empirical functions drawn on Figure 9 are simple approximations to the re- 
gression curve determined for the midpoints of the indicated range of events for 
the years 1966 through 1970. The midpoints were used as a compromise between 
the possible ways of counting multiple events. Data from 1961-1965, shown as 
open circles, were excluded from the fitting analysis because of the brief trans- 
mission schedule relative to the latter years and also the complexity in detecting 
the occurrence of events. For all years, the data were obtained by normalizing 
the observed number of events by the ratio of the total hours of recordings actu- 
ally searched (500-2000 per year) to the number of hours available for observa- 
tion (approx. 8800 and 4300 for the 24 and 12-hour functions). 

The data indicate that a 'yearly average' sunspot number near 85 the frequency of 
occurrence of SPA's over the entire globe is statistically about one per day. 

There is also the suggestion of a non-linear relationship between SPA occurrence 
and sunspot activity with the occurrence nearly vanishing during years of low 
sunspots. 

During 1961 there were 500 hours of transmission near noon while there were 
750 hrs in 1962, 700 in 1963, 600 hrs in 1964 and 1000 hrs in 1965. Thus 1961- 
1965 probabilities are less statistically certain than those based on the nearly 
continuous transmissions schedule maintained in latter years. This is especially 
true for 1961 when there were no hours of simultaneous monitoring at both ends 
of the path. It was also more difficult to measure SPA's during 1961-1965 since 
the Omega system was not operated in the modem "absolute" configuration where 
each station operates directly from cesium standards but rather in the older sys- 
tem configuration wherein some stations were operated as "masters" and others 
as "slaves", which functioned to approximate in-phase reflectors of the master 
signals delayed by retransmission location in the 10-second time-sharing commu- 
tation pattern. During much of the early period a station in the Panama Canal 
Zone was operating as master with both Hawaii and New York operating as slaves. 
SPA's for the Hawaii to New York path were derived from measurements of the 
New York - Canal Zone master-slave pair near Hawaii and of the Hawaii - Canal 
Zone pair near New York. Propagationally, the former measurement includes 
transmission from the Canal Zone to Hawaii (CH) thence from Hawaii to New 
York (HN) measured with respect to transmission from the Canal Zone to New 
York (CN): Measurement #1 represents CH + HN - CN. Similarly, Measurement 
#2 represents CN + NH - CH. Thus the sum of the two measurement represents 
HN + NH, i.e. , roundtrlp propagation between Hawaii and New York. Since 
checks during other transmission periods indicate SPA's are similar in either 
propagation direction, the derived roundtrlp is thus representative of the 
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Hawaii - New York measurement obtained directly in later years. Measurement 
of the multiple records is, however, awkward and subject to higher uncertainty 
than the direct measurement. 


Variation of SPA Magnitude with VLF Frequency 

The frequency dependence of SPA magnitudes is of interest both to the study of 
ionospheric behavior and for navigation using the Omega System. This behavior 
was investigated at the carrier frequencies of 10.2 and 13.6kHz for both the 
Hawaii to New York and Trinidad to New York paths. Data were compiled for 
approximately 50 SPA's in each of the years 1967, 1969 and 1970 for the HA-NY 
path and a total of 50 SPA's for the years 1967 and 1970 on the TD-NY path. 

The events used were selected randomly from the total set and were determined 
mainly by their availability as individual events at both carrier frequencies. 

Preliminary checks indicated that onset times and general shapes were the same 
for both frequencies for both paths being considered. The maximum phase off- 
sets then were used to determine the average and best-fit linear estimate of the 
slope of the A <j> t36 versus A <j> 1(l2 function, with the comparisons being given in 
Figure 10 and summarized in Table I. The results indicate to a high degree of 
confidence that the relative SPA magnitude effects at 13.6 and 10.2 kHz are di- 
rectly proportional, with the indicated proportionality constant of 0. 75 being 
equal to the ratio of the carrier frequency wavelengths: 

/XI 3.6 \ 

^ 13.6 “ 0- 7 5 A0 lo2 - -J Ad )0 2 (Cydcs) 

which is as expected if the phase height variations estimated from data at either 
frequency are to be the same. Since no data at other frequencies were used, the 
applicability of the wavelength-ratio proportional fy to other frequencies is 
unknown. 

Correlation of SPA Shape Characteristics 


Figures 11 through 13 show the functional interdependence of the various shape 
characteristics of the SPA's being studied. The plots are presented in a format 
intended to represent both the predictive as well as physical aspects of these 
relationships; i.e. , how well one can predict the maximum phase offset from the 
observed run-off rate, or the disturbance duration or decay rate from the ob- 
served maximum offset. These relationships have been approximated by best- 
fit (least-squares) lines through the origin and with arbitrary intercept. The re- 
sulting slopes and correlations are tabulated in Table II; ordinate intercepts are 
shown on the figures. The 'true' best-fit line was computed to show both the cor- 
rect sample correlation, r B , and the discrepancy with the line through the origin. 
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Table I 

Summary of SPA Magnitude Variation at 10. 2 and 13. 6 kHz 


PATH 

YEAR 

■ 

mSm 

A ^10.2 

AVERAGE 

AVERAGE 

SLOPE 

AVERAGE 

SLOPE 

LINEAR 

FIT 

INTER- 

CEPT 

SAMPLE 

CORRE- 

LATION 

TO— NY 

B 

47 

11.8 

B 

0.763 

0.77 

-0.1 

0.99 

HA—- NY 

1967 

48 

20.0 

14.9 

0.745 

0.65 

B 

0.96 

1969 

59 

22.5 

17.1 

0.750 

0.79 

B 

0.98 

1970 

41 

26.3 

20.4 

0.775 

0.81 

-1.0 

0.99 


AVERAGES 

0.758 

0.75 

0.0 

0.98 
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Figure 11. Maximum Phase Offset vs. Onset Phase Rate-of-Change for SPA Events 

on Hawaii to New York Path at 10. 2 kHz 
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Figure 12. Disturbance Duration vs. Maximuir 
Offset for SPA Events on Hawaii to New York Ps 
10.2 kHz 
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Figure 13. Recovery Phase Rate- or- Change vs. Maxi- 
mum Phase Offset for SPA Events on Hawaii to New York 
Path at 10.2 kHz 
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Table II 

SPA Shape Correlation 


YEAR 

NO. 

OF 

SPAs 

MAXIMUM PHASE OFFSET 
VS. 

ONSET RATE-OF-CHANGE 
• 

T o^o 

TOTAL DURATION 
VS. 

MAXIMUM PHASE OFFSET 
T ~ (l/d>)Ad> M 

RECOVERY RATE-OF-CHANGE 
VS. 

MAXIMUM PHASE OFFSET 
<A R ~0/T R )A<k 



(MIN.) 

/\ 

To 

(MIN) 

r B 

t' 

[CEC/MIN) 

/N 

, 4 > 

(CEC/MIN 

r B 

£ 

(M&.) 
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0.69 

1967 

68 

6.7 
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51.3 

0.79 

1968 

76 

5.9 
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0.53 
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46.3 

78.1 

0.65 

1969 

99 

5.4 

3.2 
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0.58 

1.33 
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82.0 

0.72 

1970 

92 

6.2 

4.7 

0.51 

0.62 

1.24 

0.65 

52.4 
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0.74 


Notes: See Figure 2 for definition of terms. Primes Indicate quantities determined by regression 
through origin (dashed lines In Figures 9-11). The quantity 4 > Is the reciprocal of the 
estimated regression slope of T on Ad* M and similarly for$ Tpand^. 




All correlation coefficients are significant in the statistical sense that the re- 
spective quantities are positively correlated at the 95% confidence level. How- 
ever, no estimated correlation slope is sufficiently high to allow explanation of 
even half of the observed scatter, although about half of the statistical variance 
can be explained in estimating recovery time. Presumably, other unobserved 
geophysical variables such as the energy spectrum of the associated flare3 or 
ionospheric chemistry variations are significant in determining the interrelation- 
ships. This statistical explanation associated with the correlation slope is in 
addition to causual explanation in the sense that when an event is recognized, the 
mean future condition is expected. Specific relations are discussed in the fol- 
lowing paragraphs. 

Maximum phase offset as a function of onset phase rate-of-change is shown in 
Figure 11. Although the indicated correlations suggest slightly shorter onset 
durations in periods of higher sunspot activity, the yearly average onset durations 
so computed yield a standard deviation of only 0. 8 cec about an overall average 
of 6. 4 cec. The apparent quantatization of rates is due to the measurement of 
onset slopes to the nearest degree. 

Total duration as a function of the maximum phase offset is shown in Figure 12. 
The relationship is primarily causual; i.e. , when an SPA has occurred, the dura- 
tion will be typically about 45 minutes (see Fig. 5). There also is a tendency for 
larger SPA’s to last longer. The median overall rate is 1.3cec/min correspond- 
ing to 0.77 minutes per centicycle of observed phase perturbation, or equivalently, 
about 6 minutes duration per kilometer of effective ionospheric height change. 
There appears to be no significant difference between years of high and low solar 
activity. 

Recovery phase rate-of-change is related to maximum phase offset in Figure 
13. The indicated correlations are higher than for the two previous functional 
relationships and are sufficient to permit explanation of half of the variance in 
recovery rates in t rms of maximum offset. The indicated recovery time is con- 
sistant with previous results (see Fig. 5) showing the range for typical-to-average 
total duration to be about 40-50 minutes. There appears to be no significant dif- 
ference b^:ween years of nigh and low solar activity. 

Spatial Variation: Normalization of SPA Magnitude for Solar Zenith Angle 


Previous sections have discussed the statistical occurence and distribution func- 
tions of SPA characteristics over a single path at one frequency and the interrela- 
tionship between SPA's at 10.2 and 13.6kHz. Spatial effects on observed SPA 
magnitude will now be studied through the use of Rome, New York phase record- 
ings of 10.2kHz transmissions from both Hawaii and Trinidad. An initial scan- 
ning indicated that the onset times and durations of events observed near noon 
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were simultaneous to within the minute or so accuracy of the instrumentation and 
that the shapes were approximately similar. Detailed attention was then directed 
to the interrelationship between peak amplitudes and illumination conditions along 
the two paths. 

Previous studies of diurnal phase variation at 10.2 kHz have shown that phase 
varies approximately proportionally to the normal component of solar flux near 
noon (although slightly modified by historical dependence). 13 Normal diurnal 
variation near noon is approximately as though the effective phase height were 
negatively proportional to the incident solar flux (cosine of the solar zenith angle) 
and the nominal phase velocity were varying negatively proportional to the phase 
height with no variation in the phase of the excitation factor associated with 
coupling-propagated energy into or from the earth- ionosphere waveguide. The 
foregoing arguments are consistent with experimental phase variation resulting 
from normal solar radiative flux variations producing ionization changes in the 
D-region near noon. We now ask whether these features are also appropriate 
for abnormal solar radiation such as produced by solar X-ray flares. The raw 
measurements were correlated to obtain the assumed linear proportionality 
between maximum phase offsets on the two paths: A <j> Trinidad - New York = 0. 44 
(A0) Hawaii -New York-0. 2 cec. The raw measurements also were normalized 
by division by the respective path average cosines of the solar zenith angles and 
essentially the same result was obtained from the correlation both numerically 
and in quality of fit. Apparently the zenith angle criterion is a poor method of 
normalization. A further comparison was made by correlating measurements 
where the average solar zenith angle over both paths was much the same so as 
to reduce complexity by minimizing differences in path illumination conditions. 

In this case, a coupling proportionality constant of 0.39 was obtained compared 
with 0.44 obtained earlier. 

The constant of proportionality between the magnitude of SPA's observed over 
the Hawaii - New York path and those observed on the Trinidad - New York path 
may be estimated theoretically. For the two paths in question secondary velocity 
variations, such as with path azimuth, are not expected to be large nor is the 
excitation factor expected to vary significantly for illumination near noon. 14 
Thus the coupling constant should be proportional tc ffective path length ratio. 
The actual total path length ratio is 0.49 = (3843 km/7814 km) but propagation 
theory developed fcr the Omega signals 4 5 recognizes a region of excitation or 
de-coupling from the earth-ionosphere wave guide of bSOkm (6.1° of arc on the 
earth’s surface) at each end. Thus the effective path length is 2 x 680km shorter 
than the geometric path length and the effective path length ratio for the two paths 
studied is (3843 - 1360 km/7814 - 1360km; - 0.39, which was the experimental 
ratio obtained when data were selected so as to minimize differences in illumina- 
tion conditions. Thus the experimental results are consistent with present theory 
on the spatial aspects of VLF propagation including the effective path length 
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concept but are not consistent with zenith angle normalization. A possible ex- 
planation may be that the phase received is not linearly related to the normal 
component of X-ray flux over a path due either to ionization equations governing 
the relationship between effective phase height and input energy or a non-linearity 
between effective phase height and velocity. In practice both may be present to 
some degree. Certainly, theory indicates marked non-linearity in the relation- 
ship between phase height and velocity for heights between 60 and 70 km, the 
range needed for the larger SPA's. Non-linearity between ionospheric phase 
height and input flux is also believed likely. Although beyond the present scope 
of effect, it would appear possible to deduce the relationship between ionospheric 
phase height and high level solar flux input through further analysis of SPA's. 
Further analysis using the methods suggested here together with the more ela- 
borate methods of Reference 4 appears particularly inviting. 


CONCLUSIONS 

The foregoing synoptic study of 500 SPA's observed between 1961 and 1970 in 
New York yields the following description of a typical SPA. A rather abrupt 
phase variation occurs reaching a maximum offset in about 6 minutes with the 
total duration of the event being about 45 minutes. Typical maximum offsets are 
about 15 cec at 10.2 kHz on the Hawaii - New York path corresponding to 15 micro- 
seconds or ionospheric phase height changes of 1. 8 km. Typical observed phase 
rates of change during onset were 3.5 cec /minute corresponding to frequency 
shift of 6 x 10' 8 or a rate-of-change of ionospheric phase height of 0.4 km/minute; 
recovery rates of change are 0.5cec/minute corresponding to a frequency varia- 
tion of 4 x 10 ' 9 and ionospheric height change rate of 0.06 km/minute. SPA's 
are observed simultaneously over the different paths and at the two different fre- 
quencies investigated. The occurrence of SPA's is related to sunspot activity 
with about one per day being observed somewhere on earth when the sunspot num- 
ber is about 85. The probability of a SPA occurring nearly vanishes during very 
low sunspot activity. Although sufficiently rare that SPA-associated phase varia- 
tions exceed nominal propagationally- induced phase scatter only about 5% of the 
time even during high sunspot activity, SPA's are of practical as well as geo- 
physical and astrophysical interest. VLF navigation, timing, and communications 
systems are expected to work safely, accurately, and reliably most of the time. 

It is primarily during unusual conditions that additional attention is needed to 
avoid false readings during these anomalies. 
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QUESTION AND ANSWER PERIOD 


QUESTION: 

I see the great correlation with sunspot activity in many other fields of science, 
but I see a weak point, and ask you what significant data do you have of any activ- 
ity that takes place on the other side of the Sun? 

MR. SWANSON: 

The solar forecasters do take into account and often will, say, call tor a higher 
probability of activity, say, 12 or 13 days after activity was last observed on the 
ground — the spot group then reappearing on the other limb. 

For the type of use that I have made of it here, I don't know if it is all that fair 
to tie it to a sunspot number directly, or whether it should be into the phase of 
a basically 11 year solar phenomenon. I wouldn't argue that. 

Really, in the short run, it is clear that if you know the types and spots and 
whether or not they are there, or whether the evolution of spots on the Sun in the 
past few weeks has gotten them to be there, you can be far more deterministic 
about it. These reports are issued regularly, incidentally, both by Boulder and 
by the Cheyenne mountain forecast. It is still, I am afraid, somewhat of a black 
art, but it is getting more refined. 

DR. WINKLER: 

It appears to me that for practical application, time service, for instance, the 
probability measure is of relatively little importance. All you can say is next 
week it is going to be bad, and you will probably have a number of disturbances, 
but you will still be unable to say when they will come and how large they will be. 

What is your opinion about comparing or correlating the various frequency dif- 
ference, or phase differences, with the view toward telling the operator at the 
moment that he has a disturbance of so many microseconds ? That i.~ really 
going back to Pierce's composite phase technique. 

MR. SWANSON: 

The correlations there are high, but not entirely perfect. I think it is worth a 
try. Let's face it, this is what we do. We shouldn't lose sight of this, it is an 
accepted technique and has been for years. If you want to adjust a frequency 
standard or something of that sort, you start making a phase track as a function 
of your local clock, and you do this more or less in the daytime. If there is a 
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sudden phase anomaly, the track makes an obvious runoff, and you say, oh, 
there was a Sudden Phase Anomaly. You wait for an hour and It goes away, and 
you are back to measuring again. 

This Is perfectly valid. Really that technique has basically the navigational 
counterpart which is to get a good dead reckoning of where you are going. We 
know from these distributions that when the event occurs the runoff would be 
fairly rapid. Your integration time is most likely on the order of a minute 
anyway. 

So, we are talking here of only a couple of integration times before we are 
beginning to get a good estimate that there has been an event. This is a very 
realistic environment. 

So, both the auto-correlated property, or sudden failure of it, is a good indica- 
tion of an SPA as well as the inner comparison in the event at several fre- 
quencies. I think most of the time we avoid these because they are so obvious. 
The danger is mostly if you are using an automatic system or something of that 
type. 


DR. REDER: 

We had made some studies of finding out whether we can predict or whether we 
can say from a measurement, let's say on NPM at 23.4 kilohertz, whether we 
can say— at Fort Monmouth— what it would be on 10.2, Haiku, Hawaii, which is 
the same path. 

And we are very happy because everything seemed to work out very fine. Unfor- 
tunately, I made some mistakes, and when I applied it to Trinidad at Fort Mon- 
mouth, it didn't work at all. 

MR. SWANSON: 

We did check the simultaneity of events on the Haiku to Rome, and also Trinidad 
to Rome path. 

DR. REDER; 

No, I am not talking about the simultaneity. This works fine. The problem was 
can it predict the size of the maximum phase, and this is the problem. 
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MR. SWANSON: 


That very subject Is in the paper. I didn't go into it here, partly because of the 
complexity. It is not as easy as you might think. It does not go strictly in a 
linear way with zenith angle, as the normal diurnal change would, near noon. 

Nor does it go strictly with the path length, apparently. 

Well, one would think the process perhaps is more easily explained than it is. 

I think there are some reasons that it doesn't work quite the way you would like 
it to. 

DR. REDER: 

Yes, but Eiic, the problem is it doesn't even work always on the same path, at 
the same time. That is the problem. 

In other words, NPM to some location and Haiku to some location. We have 
three frequencies, 10.2, 13.6, and 11-1/3 kHz. Sometimes it will work very 
well using some kind of ionospheric model, and then comes an SID and it doesn't 
work at all. 

Apparently, it means that it has a lot to do with the spectrum of the solar atmos- 
phere, and that we have to take it into account. 


MR. SWANSON: 


I am sure this is an element. I think it is probably also true that over any of 
these long paths there is quite an array of the solar zenith angles, and I think, 
as you get more of a flare flux, the depression of the ionosphere, if you will, 
and the pha?e response to this depression, both will become non-linear. 

So, now you say, I want a zenith angle, and what are we talking u'l^ut. If we are 
talking about the average, we know that is wrong all of a sudden. It is a non- 
linear effect here. I have my own speculative feelings that one might be able to 
pull inversions, and deduce the relationship in essence between say velocity or 
effective pha3e height and the zenith angle, the flux component, by running 
studies of this sort. 

In this regard, there is also the alternative approach of going directly to the 
aeronomy models using flux input on a particular flare, such at are available 
from the Solarad system. 
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OR. REDER: 


Eric, it is a great pleasure to notice that whenever you don't know an answer 
you rely on non-linearity. Welcome to the club. 

MR. SWANSON: 

I am giving two approaches. I don't know how the non-linear lies, but I do know 
there are two. If you take even a simplistic model such as Wait published years 
ago, and look at the partial of die velocity with respect to the nominal height, as 
the ionospheric height levels go down toward 60 kilometers, it becomes a much 
wider spread in velocity than is true about, say, near 80. 70 to 80 is maybe 
half the velocity change than 60 to 70 is. 

I think that the other half of the coin is what happens to the aeronomy given die 
flux coming in, and I suspect that is non-linear also. 

MR. EASTON: 

Thank you very much, Eric. 
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